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new calculations of the corresponding partonic hard-scattering functions that include the so-called 
"non-derivative" contributions not previously considered in the literature. We find a remarkably 
simple structure of the results. We also present a brief phenomenological study of the spin asymme- 
try, taking into account data from fixed-target scattering and also the latest information available 
from RHIC. We make additional predictions that may be tested experimentally at RHIC. 
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I. INTRODUCTION 

The single-transverse spin asymmetry in the process pp — > ttX is among the simplest spin 
observables in hadronic scattering. One scatters a beam of transversely polarized protons 
off unpolarized protons and measures the numbers of pions produced to either the left or 
the right of the plane spanned by the momentum and spin directions of the initial polarized 
protons. This defines a "left-right" asymmetry. Equivalently, the asymmetry may be ob- 
tained by flipping the spins of the initial polarized protons. This gives rise to the customary 
definition 



where St denotes the transverse spin vector and £ the four-momentum of the produced pion. 
We have written in short the symbol a for the cross section; we will be interested here in the 
invariant differential cross section Ed 3 a(£, st)/c1 3 £, where E is the produced pion's energy. 
We assume the pions to be produced at large transverse momentum £±. Needless to say, 
one can consider analogous single-spin asymmetries with other hadrons in the initial or final 
states, or with a final-state photon or hadronic jet. 

Measurements of single-spin asymmetries in hadronic scattering experiments over the 
past three decades have shown spectacular results. Large asymmetries of up to several tens 
of per cents were observed at forward (with respect to the polarized initial beam) angles of 
the produced pion. Until a few years ago, all these experiments were done with a polarized 
beam impeding on a fixed target (see, for example Q). These experiments necessarily had 
a relatively limited kinematic reach, in particular in £±. Now, after the advent of the first 



description is bound to be under better control. 

Despite the conceptual simplicity of A^, the theoretical analysis of single-spin asymme- 
tries in hadronic scattering is remarkably complex. The reason for this is that the asymmetry 
for a single-inclusive reaction like p^p — ► ttX (the symbol f denoting from now on the po- 
larization of the proton) is power-suppressed as l/£± in the hard scale set by the observed 
large pion transverse momentum. This is in contrast to typical double (longitudinal or 
transverse) spin asymmetries that usually scale for large £±. In essence, the leading-twist 
part cancels in the difference a(£, St) — &{£, — s*t) in the numerator of A^. That A^ must be 
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power-suppressed is easy to see: the only leading-power distribution function in the proton 
associated with transverse polarization is transversity j^|. For transversity to contribute, 
the corresponding partonic hard-scattering functions need to involve a transversely polar- 
ized quark scattering off an unpolarized one. Cross sections for such reactions vanish in 
perturbative QCD for massless quarks because they require a helicity-flip for the polarized 
quark, which the perturbative qqg vertex does not allow. In addition, a non- vanishing single- 
spin asymmetry requires the presence of a relative interaction phase between the interfering 
amplitudes for the different helicities. At leading twist this phase can only arise through a 

loop correction, which is of higher-order in the strong coupling constant and hence leads to a 

rn 

further suppression. These arguments are, in fact, more than 30 years old |6[ and led to the 
general expectation that single-spin asymmetries should be very small, in striking contrast 
with the experimental results. 

Power-suppressed contributions to hard-scattering processes are generally much harder 
to describe in QCD than leading-twist ones. In the case of the single-spin asymmetry 
in pp — > ttX, a complete and consistent framework could be developed, however p|. It 
is based on a collinear factorization theorem at non-leading twist that relates the single- 
spin cross section to convolutions of twist-three quark-gluon correlation functions for the 
polarized proton with the usual parton distributions for the unpolarized proton and the pion 
fragmentation functions, and with hard-scattering functions calculated from an interference 
of two partonic scattering amplitudes: one with a two-parton initial state and the other with 
a three-parton initial state jl, S|- As we shall review below, the necessary phases naturally 
arise in these hard-scattering functions from the interference of the two amplitudes [21, ls| ■ 
Other, related, contributions to the single-spin asymmetry have been proposed as well, 
for which the twist-three function is associated with the unpolarized proton, or with the 
fragmentation functions |9(. 

We note that also other frameworks have been considered in the literature for describing 
single-spin asymmetries in hadronic scattering. One of these introduces distribution func- 
tions that depend on intrinsic transverse momenta of partons inside the proton Q], cor- 
related with the proton spin. Because hadronic cross sections are steeply falling functions 
of £±, relatively modest intrinsic transverse momenta may generate substantial single-spin 
effects. Calculations based on this approach have had considerable phenomenological suc- 
cess [ill]; however, they rely on a factorization in terms of transverse-momentum-dependent 
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(TMD) distributions that has generally not been established so far. They should therefore 
perhaps be regarded as models for the power-suppressed An, in contrast to the framework 
developed in |7| which is derived from QCD perturbation theory 1 . 

In the present paper, we extend the work of Ref. Only a certain class of contribu- 
tions, the so-called "derivative" pieces, to be introduced in detail below, were considered 
in [3]. These indeed dominate in the kinematic regime of interest for single-spin asymme- 
tries, in particular at forward angles of the produced pion and at the lower fixed-target 
energies. Here we also derive the "non-derivative" contributions. The full structure of the 
theoretical twist-three expression for a single-spin asymmetry contains both the derivative 
and the non-derivative contributions, and it is an interesting theoretical question how the 
two contributions combine in the final result. Furthermore, from a phenomenological point 
of view, the non-derivative contributions are expected to become relevant at more central 
pion production angles and also at higher energies. 

All in all, our study is motivated to a large extent by the advent of data from the 
RHIC collider [3, 0, I^J. By establishing that large asymmetries at forward angles persist 
to high energies, measurements at RHIC have already opened a new chapter on single- 
spin asymmetries in hadronic scattering. We emphasize that at RHIC also the unpolarized 



pion production cross section has been measured, in the same 



cinematic regimes as covered 



by the measurements of the single-spin asymmetries . An overall very 

good agreement between the data and next-to-leading order (NLO) perturbative calculations 



based on collinear factorization was found 1131 1141. I15I 1. This is in contrast to the situation 
in the fixed-target regime, where a serious shortfall of NLO theory was observed ~v\ . 
Thus, it appears that the single-spin asymmetry data from RHIC, for the first time, can be 
adequately described by theoretical calculations based on collinear factorization and partonic 
hard-scattering functions calculated to low orders in perturbation theory. Even though the 
calculations described in this work are all only at the leading-order (LO) level, we are 
confident that they offer relatively solid predictions for spin asymmetries. The prospects of 
more data to come in the near-term future clearly warrant renewed and detailed theoretical 



1 We emphasize, however, that the situation is different in cases where a hard scale is present and a small 
transverse momentum is measured. Here the TMD distributions are indeed important ingredients to 
the theoretical description. For recent work on the connection between the twist-three and the TMD 
approaches for this case, see [T^ . 
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calculations and studies. We regard our paper as a significant step in that direction. 

This paper is organized as follows: in Sec. |Hl we present our calculation of the single-spin 
asymmetry in hadronic scattering. We keep the presentation as compact as possible and refer 
to Ref. 0] for some further details. We introduce the kinematics, discuss the factorization 
and then present in some detail the calculation of the partonic twist-three hard-scattering 
functions. Here we focus on the new aspect of our work, the derivation of the non-derivative 
contributions, for which we find a remarkably simple structure. In Sec. IIHI we present a 
phenomenological study using our new results. We in particular fit the unknown twist-three 
quark-gluon correlation functions to the new RHIC data and to some of the older fixed-target 
data. We use the fit results to make further predictions for spin asymmetries measurable at 
RHIC. Finally, we draw our conclusions in Sec. II VI 



II. CALCULATION OF SINGLE-SPIN ASYMMETRY 

We start by specifying our notation for the kinematics. We consider the reaction 

A(P,s T ) + B(P')^h(£)+X , (2) 

where A is a transversely polarized spin-1/2 hadron with momentum P and spin vector St, 
B is an unpolarized hadron with momentum P', and h is a hadron produced with momentum 
i. The reaction is completely inclusive otherwise. We define the Mandelstam variables 



s = 


(P- 


f P'f -2P-P' 


T = 


(P- 


- if ~ -ip ■ i , 


U = 


(P> 


- if ~ -2P' ■ e 



and the Feynman-variable 



(3) 



24 T-U 



where the last equality holds in the hadronic center-of-mass system. 



A. Factorization of the spin-dependent cross section 

As was shown in Ref. j?J], to leading power in the transverse momentum £± of the produced 
hadron, the spin-dependent cross section dAa(£_i, §t) factorizes into combinations of three- 
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field twist-3 matrix elements, twist-2 parton distributions and/or fragmentation functions, 
and partonic hard-scattering functions. The general structure of the cross section is 

A(TA+B-*hx(Z±, st) = ^i/A^i' x 2, s T ) <8> 4>b/B(x') ® H ab ^ c (£±, s T ) ® D c ^ h (z) 

abc 

+ 5^ s 9g/a(x, s t ) ® 4) ® H'ab^c{£±> s T ) ® D c ^ h (z) 

abc 

+ Yl S 9a/A(x, s T ) ® <pb/ B (x') ® H'^ c (e ± , s T ) ® D® h (z u z 2 ) 

abc 

+ higher-power corrections , (5) 

where the symbol ® denotes an appropriate convolution in partonic light-cone momentum 
fractions, to be specified below. Additional arguments, such as the pion transverse momen- 
tum or the factorization/renormalization scales, have been suppressed. The superscripts 
"(3)" in Eq. J5J) indicate the higher-twist functions. The other functions, 8(f> a /A(x) 
and D c ^h(z), are the standard twist- two unpolarized and transversity parton distributions, 
and the fragmentation functions, respectively. The sums run over all parton flavors: quarks, 
anti-quarks and gluons. As Eq. (jSJ) shows, there are in general three types of contribu- 
tions to the cross section, distinguished by the twist-3 function being associated with either 
the polarized proton (first line), the unpolarized proton (second line), or the fragmentation 
process (third line). For each of these contributions, there is a separate set of partonic 
hard-scattering cross sections, denoted by H a b_, c , H' ab _^ c , H" b _^ c in Eq. (J3J). In this paper, 
we will consider only the contributions of the first type to the spin-dependent cross section. 
The other two are expected to be suppressed relative to the first one, as discussed in |7( and 
verified by explicit calculation in jjj. 

We can think of the contribution in the first line of © in terms of the generic Feynman 
diagram shown in Fig. ^ The upper part of the diagram represents a twist-3 function 
for the polarized proton, generically given by a three-parton correlation. As was discussed 



in [7], the dominant contributions to the polarized cross section at forward production 
angles of the pion are expected from a correlation that connects two quarks and a gluon to 
the hard-scattering function. We will focus on this particular contribution as well. Other 
contributions, involving three exchanged gluons will also exist and play a possibly 
important role in production at mid-rapidity. 

In order to find the field-theoretic expression for the twist-3 function in the first line 
of Eq. ©, and to derive the rules for computing the associated hard-scattering functions, 




FIG. 1: Generic Feynman diagram contributing to the single transverse-spin asymmetry for inclu- 
sive pion production in proton-proton scattering at leading twist (twist-three). The polarized cross 
section can be factorized into convolutions of the following terms: twist-three quark-gluon correla- 
tion functions for the transversely polarized proton, parton distributions for the unpolarized proton, 
pion fragmentation functions, and hard-scattering functions calculable in Q CD perturbation theory. 

we consider the diagram in Fig. |21 Here the parts labeled T a and H a ^ c represent the 
twist-3 function and the partonic hard-scattering, respectively, which are connected by the 
two independent integrals over the momenta k\ and k 2 that they share. We thus have the 
following expression for the contribution of the diagram to the spin-dependent cross section: 



1 . f d k a Ik 

dAa(£±, st) = ^Y1 / T^i J^i T a( k i, k 2 , s T ) Hab^ c (ki, k 2 , £±) <g> 0&/b(z') ® D c ^ h {z) , 

abc 

(6) 

where 1/2S is a flux factor and the sum again runs over flavors. In the above expression, 
spinor, color and Lorentz indices connecting the hard and long-distance parts have already 
been separated, using the techniques developed in 0, as sketched in Fig. El In a covariant 
gauge, the function H ab ^ c (ki, k 2 ) is contracted with ^ N -t_i ) (t B )ij[(^ ItypP^I '(27r) , where 
the factor (2tt) is due to the normalization of the twist-3 matrix element T a , Nc = 3 is the 
number of colors, B and i,j are the color indices of the initial gluon and quarks, respectively, 
and the matrices {t B )ij are the SU(3) generators in the fundamental representation. 

The next step is to perform a "collinear" expansion of the expression for the diagram . 
Due to perturbative pinch singularities of the partonic scattering diagrams jj, the integra- 
tion in Eq. (jUJ) is dominated by the phase space where kf ~ 0, and we can approximate the 
parton momenta hi entering the hard scattering to be on-shell and nearly parallel to the 
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FIG. 2: Factorization of hard part and twist-three matrix element: (a) before and (b) after sepa- 
ration of spinor trace and Lorentz indices. For simplicity we have omitted the unpolarized parton 
distribution and the fragmentation function, n ~ P' is a light-like vector with spatial components 
in direction opposite to those of the initial momentum P. 



momentum P of the initial polarized proton: 



kf = Xi P» + k» + kl T /( Xl S)P'» 



(7) 



where the are perpendicular to both P and P' , and where kf T = — (/cf ± ) 2 . The last 
term oc k? ± in ((7j) can be neglected since it is beyond the order in fa t ± that we consider. 
The collinear expansion enables us to reduce the four- dimensional integrals in Eq. © to 
convolutions in the light-cone momentum fractions of the initial partons. Expanding H a b^, c 
in the partonic momenta, fa and k2, around k\ = x\P and fa = x 2 P, respectively, we have 



H ab ^ c (fa, fa) = H a 



,(xi,x 2 ) + d^ b ~" c (xi,x 2 ) (fa -xxPf 
oki 

+ dH ^7 C (xi, x 2 ) (fa - x 2 Pf + ... . 
ok 2 



(8) 



Because of (|7jl. the derivatives in the latter equation are in the transverse vectors k^ ± only. 
The expansion (|SJ), substituted into Eq. ©, allows us to integrate over three of the four 
components of each of the loop momenta fa. The top part of the diagram T a then becomes 
a twist-three light cone matrix element, given by 



T a Axi,x 2 ) 



dU\ dy 2 c i Xl P+y- +i(x 2 -x 1 )P+y- 
47T 

x (p, s T fe(oh + b STann F a + (y 2 )] 4>a(yT)\P, 5r> 



(9) 



where we have introduced the subscript U F" to indicate that the matrix element involves the 
gluon field strength tensor F a + . The additional ordered exponentials of the gauge field that 
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(c) 

FIG. 3: Three classes of quark-gluon scattering diagrams contributing to the spin- dependent cross 
section A<t(st): (a) diagrams with an initial-state pole, (b) and (c) diagrams with a final-state pole. 
Symbols B and ij are color indices for the gluon and the quarks. The propagator that provides the 
pole is indicated by a bar. All poles shown are "soft-gluon" poles, contributing at x\ = x 2 (see text). 

make this matrix element gauge invariant have been suppressed Eq. (jHJ) as it stands 
is valid in the light-cone gauge A + = 0. T a ^p is a symmetric function of its arguments, 
T atF (x 1 ,x 2 ) =T a>F (x 2 ,x 1 ). 



B. Poles in hard-scattering functions and contributions to fc^-expansion 

In addition, T a> p is real, implying that the phase needed to generate a single-spin asym- 
metry has to arise in the functions H a ^ c in Eq. (jHj). As was shown in [?], imaginary 
parts in H a b^ c can arise even at tree level, thanks to the pole structure of the hard scatter- 
ing function. Examples of this are shown in Fig. El for the case of quark-gluon scattering. 
Imaginary parts arise from the scattering amplitude with an extra initial-state gluon when 
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its momentum integral is evaluated by the residues of unpinched poles of the propagators 
indicated by the bars. The on-shell condition associated with any such pole fixes the mo- 
mentum fraction of the extra initial-state gluon and hence further simplifies the integrations 
over the momenta k\ and k 2 in Eq. (jUJ). Roughly speaking, all of the diagrams in Fig. |3] 
provide an unpinched pole at X\ = x 2 , with subtleties that we will address shortly. At these 
poles, one has [tJ 

dH ab ^ c dH ab ^ c 
dk p (x 1 ,x 2 = x 1 ) = - (x 1 ,x 2 = x 1 ). (10) 

Thanks to this property, one can organize the calculation of the partonic hard-scattering 
functions with a simpler momentum flow, using a single transverse momentum k±, as shown 
in Fig. |H 

In order to demonstrate the emergence of a strong-interaction phase through a pole 
contribution at x± = x 2 , let us consider the specific example for the initial-state interaction 
shown in Fig. Ufa). We need to consider contributions for which the initial-state gluon 
attaches on the right or on the left side of the cut. The propagator denoted by a bar in the 
left part of the figure reads 

1 1 



(x'P' + (x 2 — X\)P + k±) 2 + ie (x 2 — Xi)x'S + ie 



0(k 



1 7T 

^-^S(x 2 - Xl ), (11) 

where in the second line we have extracted the imaginary part provided by the propagator, 
which contributes to the single-spin asymmetry. When the gluon attaches on the right-hand 
side of the cut, we obtain the same result, but with opposite sign. Therefore, effectively the 
difference of the two diagrams in Fig. Ufa) contributes. The first term on the right-hand 
side of Eq. (JHJ) cancels in this difference, so that eventually only the other two contribute. 

In Eq. pijl we have neglected a term cx k\ since we are only interested in first-order 
(linear) k? effects. A linear term is not present in the delta-function in (fTT|) because the 
vector k± is perpendicular to both P and P' . For final-state interaction, this situation 
changes. Generic diagrams with final-state interactions involving the "observed" parton are 
shown in Fig. Efb). On the left side of the diagram, a phase from the propagator marked 
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(6) 



FIG. 4: Specific examples of diagrams for generic (a) initial-state and (b) final-state interactions, 
along with simplified notation of external momenta. 



by the bar arises as 



1 



(£/z - (x 2 - Xi)P - k ± ) 2 + ie -2P ■ £ (x 2 - x x )/z - 2£ ■ k ± /z + ie 



+ 0{k 2 T ) 



mz 



p 6{x 2 -x 1 -2£-k ± /T) , (12) 

where the momentum of the fragmenting ( "observed" ) final-state parton is related to that of 
the produced hadron by £ = zp c , and T has been defined in Eq. (J3J). Again, the propagator 
on the right side of the cut has the same pole, with opposite sign. As Eq. (fT2"|) shows, the pole 
provided by the final-state interactions is located near x\ = x 2 , but displaced by a term linear 
in k±. When inserted into the collinear expansion JHJ), this term will make a contribution to 
the single-spin asymmetry involving a derivative of the delta-function and hence, by partial 
integration, a derivative of the twist-three quark-gluon correlation function 

Such "derivative" terms may, however, also arise in a different way, through the on-shell 
condition for the unobserved final-state parton. For the diagrams on the left-hand-side of 
Figs. Efa) and (b), the momentum carried by that parton is p% = x'P' + x 2 P + k± — £/z, 
where the superscript L (R introduced later) refers to the diagrams whose extra initial-state 
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gluon is attached to the left (right) of the cut. The phase space provides a delta-function 
that puts this particle on its mass-shell, 

5{(p L d f) = 5(x 2 {x'S + T/z)+x'U/z-2£-kJz) 

6 x 2 - x ~ . c I rr, > (13) 



x'S + T/z \ x'S + T/z 

where 

- X ' U I Z (<iA\ 

x = „ . 14 

x'b + T/z 

x can be interpreted as the "usual" value of the partonic momentum fraction of the polarized 
proton if there is no k±. Eq. (|T3*j) fixes x 2 in terms of the Mandelstam variables S, T, U and 
a linear term in k±. The latter will give rise to "derivative" contributions in the same way 
as Eq. (jTJJ) does. If the initial gluon attaches on the right-hand-side of the cut, however, the 
momentum of the unobserved parton is = x'P' + x\P — t/z, and the resulting on-shell 
condition fixes x\\ 

S ((pf ) 2 ) = 5( Xl - x) , (15) 

with no dependence on k±. 

Additional contributions to the collinear expansion can of course also arise from terms lin- 
ear in k± in the other "hard" propagators or in the numerator of each diagram. These terms 
do not lead to "derivative" contributions to the single-spin asymmetry, but to contributions 
involving T a ^p itself. 

We close this section with two further observations. First, we note that final-state inter- 
actions involving the "unobserved" parton d cancel when summing over contributions where 
the additional gluon attaches on the right or the left side of the cut. Second, the contribu- 
tions we have discussed are all characterized by the additional initial gluon becoming soft. 
The poles arising from this are, therefore, customarily referred to as "soft-gluon poles" . The 
hard-scattering diagrams will in general possess also other poles, for which an initial quark 
becomes soft |7(. Such "soft-fermion poles" are expected to play a less important role and 
are not considered in this work. 

In the next section, we will provide a "master formula" that allows to take into account 
all contributions to the /cj_-expansion discussed above simultaneously and in a systematic 
and relatively straightforward manner. 
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C. "Master formula" 

As was shown in jjj, the factorized expression for Acr(sr) takes the form 

dAa(s T ) oc -^^2 J dzD c ^ h (z) J ^j- <p b/B (x') J dx 1 dx 2 T ajF (x 1 ,x 2 ) 

abc 

xie pSTnn hm ■^ F H^ c (x l ,x 2 , x' , z) , (16) 

where 

e ps T nn = e p*nu ST ^ n ^ 

with n and n two light-like vectors whose spatial components are parallel to those of P' and 
P, respectively. According to the discussion in the previous section, we are therefore led to 
consider the following general expression: 



lim -^-p I dxx I dx 2 T a , F (x u x 2 ) 
k ± ^o ok. 



H L (xi,x 2l k±) 8(xi — X2+VX- k ± )5(x 2 - x - v 2 ■ k±) 
—H R (x 1 ,x 2 ,k_ L )5(x 1 —x 2 +Vi-k±)5(xi—x) . (18) 

Here, Hl and Hr denote the contributions to H a ^ c for any diagram, when the initial gluon 
attaches on the left or the right side, respectively. v% and v 2 are vectors made of P, P', 
and i whose form follows directly from the preceding discussion. The first delta-function in 
each of the two terms associated with Pl,r results from the propagator poles discussed in 
Eqs. (fTTj) and (fI2j). For initial-state interactions, v\ — (see Eq. (fTTjl). for final-state ones, 
V\ = 2£/T (see Eq. ()12)) ). An important point is that v% is the same vector on the left and 
on the right-hand-side of the cut. The second delta-function in each term results from the 
on-shell condition for the unobserved particle. As explained earlier, these delta functions 
differ for the two sides of the cut. We have v 2 = 2£/(x'zS + T) for the left side and v 2 = 
for the right one, which we have already used. 

A straightforward way of dealing with the expression in (JTHj) is to use the various delta- 
functions to perform the integrations over X\ and x 2 . This gives: 





lim 

fe±^o dk 



p 



T a ,F(x + (t> 2 - vi) ■ k±, x + v 2 ■ k±) H L (x + (v 2 - v x ) ■ k ± , x + v 2 ■ k ± , k±) 
-T a>F (x,x + vi ■ k±) H R (x,x + v 1 ■ k±, k±) . (19) 
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This term can be organized as 

(v 2 ~ v 1 ) p H L (x, x, 0) dTa ^ X) + T a>F (x, x) x (20) 



r 9 
x nm 



fcx^o dk p . 



J fci =o 



H L (x + (v 2 - v\) ■ k±, x + v 2 • k±, k±) - H R (x } x + vi ■ k±, fcj 

This is our "master formula" . In deriving it we have used that the hard-scattering functions 
with the gluon attaching on the left or the right side of the cut are the same at k± = 0, 

H L {x,x,0) = H R {x,x,0) . (21) 

Equation (|2()|) applies to both initial- and final-state interactions. As one can see, the 
first term is proportional to dT a ^(x,x)/dx [thanks to Eq. ()21|) it does not matter whether 
we write Hl(x,x,0) or H R (x,x,0) in this term]. This is the "derivative" contribution that 
we discussed above and that was originally computed in Ref. The second term involves 
only T ai p(x, x), without a derivative. Equation (J2U|) allows a simultaneous computation of 
both the derivative and non-derivative contributions. 

The next step is to consider all contributing partonic reactions and to calculate the contri- 
butions to Eq. (J20J). The partonic channels we need to consider are (qg)g — > qg, (qg)g — > gq, 

(qg)q 99, (qgW -> w', (qg)q' -> q'q, (qg)q -> qq, (qg)q -> q'q', (qg)q -> ?g, where for 

each the first two initial partons are entering from the polarized proton via the twist-three 
correlation function T a> p. We remind the reader that we are ignoring contributions involving 
a three-gluon twist-three correlation function, which would correspond to a (gg) initial state. 
Crossed channels are implicit and taken into account as well. 

Upon calculating all associated hard-scattering functions, we found that they possess a 
remarkable property: for each process, 

d r 1 

H L (x + (v 2 - vi) ■ k±, x + v 2 ■ k±, k±) - H R (x, x + vi ■ k±, k±) 

J k ± =0 

H L (x,x,0), (22) 

which is again valid for the case of both initial- and the final-state interactions. We have 
not been able to develop a proof why Eq. (j22j) holds in general, even though the equation 
is certainly not accidental and such a proof should be possible. In any case, Equation ()22|) 
leads to a dramatic simplification of the final result. Inserting (J22*|) into (j20j) . one finds the 
expression 

( '" )p H L (x,x,0) [T a , F (x,x)-xT^ F (x,x)] , (23) 



lim 

fcx^O ok 



p 



(V2 - V X )p 



X 
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where we have used the short-hand notation T^ F (x, x) = dT at p(x,x)/dx. Thus, even 
though there could have in principle been two separate hard-scattering functions multiply- 
ing T a ^p(x,x) and T' aF {x, x), the final result for the combined derivative and non-derivative 
terms will have a single hard-scattering function for each process, summed over initial- and 
final-state contributions and multiplying simply the combination T ai p(x,x) — xT' aF {x : x). 
This hard-scattering function is furthermore identical to the one calculated for the deriva- 
tive p.ece in Q. The ™ structure is tueu very akl u to that of the unpolled cross 
section. We are now in the position to give the final answer for the single-spin asymmetry. 



D. Final result 



For definiteness, we recall the expressions for the vectors v i and v 2 introduced above. We 
have 

21 2 Pc x 

U2 = ^sTt = -^' (24) 

with the partonic Mandelstam variable u = (p c —p') 2 = x'U/z. For initial-state interactions, 
see Eq. (fTTj) . we have V\ = 0, while for final-state ones, see Eq. (|T2j) . 



Vl= T 



21 2p c x 



t 



or t> 2 — vi 



2p c x 



u 



t 



(25) 



where i = (p c — p) 2 = xT/z. 



Using (|T6*j) and following the steps presented in detail in Ref. jjj, we then find the final 
expression for the polarized cross section: 



d 3 A. 



a{s T ) 



dH 



a t f dz f dx' 1 , 

a,b,C Zraln X min 



(26) 



T a p(x, x) - x\ —T„ F (x, x 



d 



zu J x 1 \dx 

where x has been defined in Eq. (JHJ), and where 

, -T/z T + U 

Xmin ~ S + U/z ' Zmin ~ S~ 

The Hab^c are the final hard-scattering functions and read 



H ab ^ c (s,i,u) , 



(27) 



H, 



ab- 



H^ c (s, t, u) + H^ c (s, t, u) ( 1 + 7 



(28) 
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where H I ab ^ rC (H^ b ^ c ) denote the contributions due to the initial-state (final-state) interac- 
tions. The factor (1 + u/t) results from the expression for t>2 — V\ in Eq. ()25|l. We have 
collected all H^ b _^ c and H^_^ c in Appendix A. Thanks to the structure we have found, the 
must coincide with the hard-scattering functions calculated for the derivative part in Ref. 
which they do, up to trivial corrections we found for some of the color factors in The 
results presented in Appendix A are also in a more compact and transparent notation. 

We emphasize again the simplicity of the structure in Eq. (|26p. which is very similar to 
that of the unpolarized cross section in the denominator of the spin asymmetry. The latter 
reads: 

(f'o y-> r 1 dz _ /'' d.v' 1 

a,b,c Zmin 



,_, d 3 a a? v-^ f 1 dz ^ . , f 1 dx' 1 . Is 

% = iEl.^"*)| If IFsTtJI 1 (29) 



1 



x -cj) a / A {x)H^ b _ c {s,t,u) 



X 



with unpolarized hard-scattering functions H^ b _^ c and the usual unpolarized parton distribu- 
tion functions in hadron A, (j) a /A{x). We give the well-known 20] H^ b _ iC also in Appendix A. 

We finally note that we have written the hard-scattering functions for both the spin- 
dependent and for the unpolarized case as dimensionless functions. The power-suppression 
of the single-spin asymmetry is then explicitly visible by the denominator u in Eq. ()26|) . Fur- 
thermore, note the factor ^/q~ s in that equation, which results from the additional interaction 
with a gluon field in the hard-scattering functions for the single-spin case. 

III. PHENOMENOLOGICAL STUDY 

We now present some first numerical results for the single-spin asymmetry derived in the 
previous section. We do not aim at a full-fledged analysis of all the hadronic single-spin 
data at this point, but would like to examine a few of the salient features of the new RHIC 
data and of the earlier E704 fixed-target pion production data. We reserve a more detailed 
analysis to a future publication. 

Let us begin by specifying the main ingredients to our calculations. We first remind the 
reader that all our calculations of the hard-scattering functions are only at the LO level. We 
therefore use LO parton distribution and fragmentation functions throughout, as well as the 
one-loop expression for the strong coupling constant. For the unpolarized cross section we 



use the LO CTEQ5L parton distribution functions 
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2l| . Our choice for the fragmentation 



functions are the LO functions presented in Ref. 22]. These have the advantage that they 
provide separate sets for positively and negatively charged pions, which are needed for the 
comparison to the experimental data. 

For the present study, we will make rather simple models for the twist-three quark- 
gluon correlation functions T a ^p(x, x) (a = u, u, d, d, s, s), relating them to their unpolarized 
leading-twist counterparts. We recall that we do not include any purely gluonic twist- 
three correlation functions, even though we do take into account the gluon-gluon scattering 
contribution in the unpolarized cross section in the denominator of A^. Our ansatz for the 
correlation functions is simply: 

T a>F (x, x } fi) = N a x a « (1 - xf« <f> a {x, y) , (30) 

where <j) a (x,fi) is the usual twist-two parton distribution for flavor of type a in a proton. 
Note that we have now written out the dependence of the functions on a factorization scale /i, 
which we will always choose as \x = £±. We will in fact assume that the functions T a} p{x, x, /i) 
evolve in the same way as the corresponding unpolarized leading-twist distributions. This 
will certainly not be correct in general, because of the different twist of the two types of 
distributions, but may be hoped to be a reasonable assumption at the moderate to relatively 
large x we are interested in here. 

We determine the parameters in Eq. (jHOj) through a "global" fit to experimental data for 
A N as functions of xp defined in Eq. (HJ), using the expressions in Eqs. (J2fijl and (J2HJ). Here 
we choose the fixed-target scattering data at \J~S ~ 20 GeV by the E704 experiment ^| for 
p^p — > ti 0,± X and p^p — > 7r 0,± X, and the latest preliminary RHIC data at y/S = 200 GeV 
by the STAR [J [for p^p -> vr^X] and BRAHMS [for p^p -> n ± X and p^p -> K±X] 
collaborations. The perturbative hard-scattering expression we have derived in the previous 
section is expected to be only applicable at high transverse momentum, starting from £± > 
a few GeV. Most of the available data points for An are at ia-values not much greater 
than 1 GeV, however. In case of the RHIC data, we always use the correct value for £± 
for each data point, keeping however only points with £j_ > 1 GeV. For the E704 data the 
situation is more complicated as most of the data points have £± ^$ 1 GeV. In addition, as 
we discussed in the Introduction, there is generally a problem with the description of even 
the unpolarized cross sections in the fixed-target regime, when hard-scattering calculations 
at low orders of perturbation theory are used. All-order resummations 



171 ] may be very 
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relevant here, which are likely to affect the spin-dependent and the unpolarized cross section 
in different ways. In view of this, we are tempted to exclude the E704 data from our analysis. 
On the other hand, the information on single-spin asymmetries is overall still rather sparse, 
and any information is potentially helpful. In particular, data for anti-proton scattering are 
only provided by the E704 experiment. Therefore, in order to include the E704 data in the 
fit, we choose £± = 1.2 GeV for these data. In addition, we allow a large shift of the overall 
normalization of the theory result used for the comparison to these data. This shift is meant 
to represent in particular the possibly large higher-order effects on A N just described. 

We have performed two separate fits to the data. One is a "two- flavor" fit, for which we use 
only the two valence densities u v and d v in the ansatz (|30|) and set all other distributions to 
zero. For this fit we introduce a normalization factor N-^ mi = 0.5 for the theory asymmetries 
in the kinematic region of the E704 data and find: 

Fit I : N Uv = 0.275 , N dv = -0.365 , 

a Uv = 0.508, (3 Uv = 0.399, a dv = -0.108 , (3 dv = 0.287 . (31) 

The fit has a x 2 -value of 304.6 for the 60 data points and is therefore of rather poor quality. 
Nonetheless, as one can see from the comparison of the fit to the experimental data shown 
by the solid lines in Fig. EJandlHl it has overall the right qualitative features. For example, 
for p^p — > tt + X the asymmetry is positive, which is reflected in a positive valence-w twist- 
three correlation function emerging from the fit. Likewise, the fact that is negative 
for p^p — > 7i~X implies a negative valence-ci distribution. For anti-proton scattering, the 
respective asymmetries are then necessarily opposite, because one has Q| 

yanti— proton ^-iproton (32) 

The asymmetries for ir° production are between those for tc + and it~ . The same qualitative 
features persist to RHIC energies, as can be seen from the comparison to the STAR (tt°) 
and BRAHMS (vr ± ) data in Fig. El 

For the second fit, we allow also sea- and anti-quark Tp functions. We then find the 
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following parameters: 



Fit II : N Uv = 0.353 , N dv = -0.594 , 

N u = -N S = -N u _ = N d _ = -19.8 , N S = N- S = -6.63 , 
a Uv = 0.696, (3 Uv = 0.559, ^ = 0.312, (3 dv = 0.488 , 

Ofi = ««sca = «<? = «4ca = «S = «S = 2.91 , 

Pu = Pu sc , = P d = Pd SEa = Ps = P* = 0-351 . (33) 

Here the relations among the various parameters for sea and anti-quarks are not fit results, 
but have been imposed. As before, we have a normalization factor A^e 70 4 = 0.5 for the 
calculated theory asymmetries at the E704 kinematics. The results of this fit are also shown 
in Figs. El and by the dashed lines. One can see that the fit is rather similar to Fit I, 
but does slightly better. Indeed, the fit has x 2 = 292.6. While the valence-quark densities 
completely dominate for the fixed-target case, the sea distributions play a somewhat more 
significant role at RHIC. We note, however, that in the case of An in K~ production (see 




FIG. 5: Comparison of the single-spin asymmetries An using our fit results in Eqs. 
to the data from E704 £3/. The solid lines are for Fit I (Eq. and the dashed ones are for 

Fit II (Eq. The lower dotted lines in the upper left part of the figure show the contributions 

to An for 7r ± production by the "non- derivative" terms alone, for Fit I. Note that the theory curves 
in the figure are normalized by N^yoa = 0.5. 
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FIG. 6: Comparison of the single-spin asymmetries An using our fit results in Eqs. \3b\) . Hn\) . \33\) 
to the preliminary RHIC data by the STAR j^/ (left) and BRAHMS jj/ (right) collaborations. The 
solid lines are for Fit I (Eq. I.VJj) ). and the dashed ones are for Fit II (Eq. 4.V.V)) ). The lower dotted 
line in the figure shows the contribution to An by the "non- derivative" terms alone, for Fit I. 

Fig. EJ) even our fit with a sea distribution does not lead to a significant change in the 
theoretical result. This is surprising at first sight, because the K~ has no valence quarks 
in common with the proton, so that sea quarks and anti-quarks should be particularly 
important here. We found that the precise admixture of valence ("favored"), non- valence 
( "un- favored" ) , and gluon fragmentation functions is very relevant in this well as that 

of the hard-scattering functions. We could improve the description of An in K~ production 
only by assuming a very large negative correlation function T^,f- 

We also address the numerical relevance of the "non-derivative" terms that we have 
calculated in this work. The dotted lines in the upper left part of Fig. El and in the left part 
of Fig. El show the contributions to An that one obtains from the "non- derivative" terms 
alone, for the case of the two-flavor fit (Fit I). One can see that these contributions are of 
relatively moderate (~ few %) size, but non- negligible. They play a bigger role at RHIC 
energies. There is roughly a 25% increase in the value of \ 2 when the "non-derivative" 
contributions are neglected. Of course, one could refit the T a> p distributions without the 
"non-derivative" contributions, in which case the theoretical spin asymmetry would be again 
very close to the dashed or solid lines in the figures. However, we found that such a fit has 
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FIG. 7: T a ^p distributions for a = u v ,d v ,u,d resulting from our fits in Eqs. and at scale 

[i = 2 GeV. We also show the corresponding unpolarized parton distribution functions )2&} . scaled 
by 1/10. 

a slightly worse y 2 , and in any case it leads to a fairly different set of T a p distributions. 

In Fig.[7|we show the T a> p distributions that we have found in our fits, for u and d valence- 
and anti-quarks, at scale fi = 2 GeV. The dashed lines are for the "two-flavor" Fit I, while 
the dotted ones are for Fit II. As one can see, the valence distributions are rather similar 
in the two fits. Only Fit II has anti-quark distributions. For all distributions, we also show 
the corresponding unpolarized leading-twist densities, scaled by 1/10 for better visibility. 

It is interesting to speculate about the reasons why the overall quality of our fits is 
relatively poor. We first remind the reader that for the reasons discussed earlier we have 
rescaled all theory asymmetries in the kinematic region of the E704 data by a factor 1/2 in 
the fit. Without the rescaling factor, the total x 2 of the fit would be increased by almost 100 
units from the current ~ 300, while the sign and the general shape of the asymmetries would 
still be consistent with the data. Small changes in the normalization of the RHIC data sets 
do not lead to a very significant further reduction of \ 2 ■ We also found that an even better 
description of all RHIC data is possible if one excludes the E704 data from the fit. Such a 
fit then tends to badly describe the data from E704, even when a normalization factor 
is applied to the latter. We recall once more that the E704 data are in a kinematic regime 
where the theoretical calculation of even the unpolarized cross section is challenging, and 
that we set £± = 1.2 GeV for them. It is therefore perhaps not surprising that we find that 
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the consistency of the total data set for An appears limited. We also caution the reader, 
however, that much of the RHIC data is still preliminary and one needs to await further 
experimental information before drawing final conclusions. 

We now use our fitted twist-three correlation functions T a ^p of Eqs. (j30j) . (|31j) . and (J33j) 
to make a set of further predictions that may be tested at RHIC. The first one concerns the 
dependence of An on the produced hadron's transverse momentum £±. This is a particularly 
interesting observable, given the power-suppressed nature of A n- In fact, as we discussed 
in the Introduction, A N is expected to decrease as l/£±, at a given xp. In Fig. |H^a) we 
plot An for tt° production at = 200 GeV at three fixed values of the Feynman variable, 
xp = 0.35, 0.45, 0.55, for our two sets of T a p in Fit I and Fit II. One can clearly see the 
fall-off with £j_. In order to experimentally verify this fall-off, that is, to keep xp fixed while 
varying £±, one would need to vary the scattering angle. On the other hand, if measurements 
are made at a (roughly) fixed scattering angle 9 or pseudo-rapidity r\ = — lntan(6 l /2), xp 
will increase along with £j_, as seen from the relation 

2£ ± 

x F = -j= sinh(?7) . (34) 

This is often the experimentally more relevant situation. As one can see from Fig.|Sfa), even 
though An decreases with £± at fixed xf, its increase with increasing i^ata given £± appears 
to be stronger. Hence, one expects that for measurements at fixed scattering angle An will 




FIG. 8: (a) Dependence of An for ir° production at RHIC at \/~S = 200 GeV on £±, for three 
different values of xp. Solid lines are for the T aj p distributions of Fit I, dashed ones are for Fit II. 
(b) Dependence on £± for fixed pion pseudo-rapidity n = 3.8 and when taking an average over the 
bin 3 < rj < 4. 



22 




60 

E„/GeV 

FIG. 9: (a) Dependence of An for ir° production at RHIC at = 200 GeV on the pion energy 
E-w, for fixed r/ = 3.3 and rj = 3.8. Solid lines are for the T a p distributions of Fit I, dashed ones 
are for Fit II. 

actually increase with £±. Indeed, this is the shown in Fig. IBfb). If one averages 

experimentally over a bin of forward rapidities, say, 3 < rj < 4, the increase of with £± 
is less pro nounced but still there. Another observable, most relevant to measurements at 



STAR |13|, is the dependence on the asymmetry on the pion energy E n = £± cosh(?y) for 
fixed rj. This is plotted in Fig. |U]for our two fits, for the cases rj = 3.3 and rj = 3.8. 

It is also interesting to consider the energy dependence of A N . So far, we have data 
from fixed-target scattering at y/S ~ 20 GeV and from RHIC at a center-of-mass energy 
about an order of magnitude higher. In order to shed further light on the mechanisms 
responsible for the large observed values of A^ at these two energies, information at an 
intermediate energy will be particularly useful. In the 2006 run, data have been taken at 
RHIC at V5 = 62.4 GeV. Using our above fit results, we find the theoretical expectations 
for A^ for 7r =l= and 7r° production as functions of xp shown in Fig. We have for now again 
correlated xf and £± through Eq. (p?4^ at fixed pseudo-rapidity rj = 3.3. Clearly, comparison 
to eventual data will require implementation of the correct kinematics. In the figure, we 
compare results for v^S* = 62.4 GeV and 200 GeV. One can see that, at fixed Xp and rj, a 
significant increase of A^ with energy should be expected. 

We finally briefly turn to processes other than inclusive-hadron production. We first 
consider the spin asymmetry in single-inclusive jet production. The partonic hard-scattering 
functions for this case are the same as for hadron production, but there are no fragmentation 
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FIG. 10: Comparison of An for ir^ and tt° production as functions of xf at y/S = 62.4 and at 
200 GeV, using our two fit results. We have chosen a fixed pion pseudo-rapidity r/ = 3.3. 

functions here. The result for An for jet-production at RHIC at forward Xp and fixed 
7] = 3.3 is shown in Fig. ^2 For comparison, we also show again the corresponding curve 
for 7T° production for Fit I. One can see that essentially the asymmetry for jets is shifted 
by a factor ~ 2 to the right with respect to that for pions. This can be understood from 
the fact that in the kinematic regime relevant here a pion takes on average roughly 50% of 
a fragmenting parton's energy [23], whereas all of the energy goes into a jet. 

Another process of interest is prompt-photon production. Photons are generally much 
less copiously produced at RHIC energies than pions, which results in larger statistical 
uncertainties on the spin asymmetries. However, given the progress on luminosity and po- 
larization at RHIC, first measurements of An for prompt photons should become possible 
in the near-term future. Photons have the advantage that one important production mech- 
anism is quark-gluon Compton scattering, qg 79, with the reaction qq — > 7$ yielding a 
smaller contribution. In addition, in these processes the photon couples in a "direct" (or, 
point-like) way, that is, there are no fragmentation functions involved. Photons can, how- 
ever, also be produced in jet fragmentation j^. The relative importance of the "direct" 
and the fragmentation contributions depends on kinematics, but also on aspects of the ex- 
perimental measurement. It is possible, for example, to largely suppress the fragmentation 
contribution by a so-called photon isolation cut j24J . In the following, in order to obtain first 
estimates, we will calculate the single-spin asymmetry for prompt photons based on either 



24 



0.08 

A N 

0.06 
0.04 
0.02 
0.00 

-0.02 

0.2 0.4 0.6 0.8 

x F 

FIG. 11: Single-spin asymmetry for jet production (with transverse momentum l± > 3 GeV) at 
RHIC at VS = 200 GeV, as a function of xp for fixed pseudo-rapidity n = 3.3. We show the results 
for both Fit I (solid) and Fit II (dashed). The dotted curve shows the same for ir° production for 
Fit I. 

the "direct" contributions alone, or on the sum of the "direct" and the full fragmentation 
contributions. The former is more representative of the asymmetry for an isolated photon 
cross section, while the latter corresponds to a fully inclusive measurement. When data will 
become available, a more careful theoretical analysis will clearly become necessary. 

Predictions for for prompt-photon production can then be obtained from Eq. (j2T?j) by 
using the appropriate hard-scattering functions for the reactions (qg)g — > jq and (qg)q — > jg 
and by replacing the coupling factor a 2 by a s a e , m , e 2 , where a e .m. is the electromagnetic 
coupling constant and e q is the fractional electric charge carried by the quark of flavor q. 
We give the resulting expressions in Appendix B, along with the corresponding ones for the 
unpolarized case, to be used in Eq. with the same replacement of the couplings. Taking 
the fit results of Eqs. (|30 p .(|5I )) .(|33 )l we obtain the predictions shown in Fig.[T2l which are for 
\/~S = 200 GeV. Again we have chosen a fixed pseudo-rapidity r\ = 3.3. For comparison we 
also show again the corresponding result for An for 7r° production. The differences between 
the asymmetries for photons ("direct only") and 7r° are quite striking. They mostly result 
from a rather different structure of the corresponding hard-scattering functions (see Eqs. (|3T)|) 
and in Appendices A and B, respectively) and are therefore a real prediction of the 
formalism. One also sees that the two fits I and II give somewhat different predictions for 
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An- for the photons "direct only" case. This is due to the contributions from (qg)q scattering 
which are present for Fit II, but absent for fit I since for this fit we assumed that the sea 
quark T a> p functions vanish. 

When the fragmentation contribution to the prompt photon cross section is taken into 
account, the single-spin asymmetry becomes much more like the one for tt°, at least for 
the lower xp. The reason is that in the kinematic regime relevant here, i.e. at relatively 
low transverse momenta £±, the fragmentation component actually dominates the cross 
section. Future measurements of the single-spin asymmetry for isolated photons should see 
an asymmetry close to the lower one ( "direct only" ) in Fig. El while for the fully inclusive 
(non-isolated) case A N should be smaller and closer to that for tt° production. 

IV. CONCLUSIONS AND OUTLOOK 

We have presented a new study of the single-spin asymmetry in single-inclusive hadron 
production in hadronic scattering. The importance of this asymmetry lies in the new insights 
into nucleon structure it may provide, but also in the challenge that its description poses for 
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FIG. 12: Single-spin asymmetry for prompt-photon production at RHIC at y/~S = 200 GeV, as a 
function of xf for fixed pseudo-rapidity r/ = 3.3. We show the predictions for both Fit I (solid) and 
Fit II (dashed). We show separately the results for the cases when the fragmentation component 
is taken into account or neglected. The dotted curve shows the earlier result for tt° production for 
Fit I. 
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QCD theory due its power-suppressed nature. We have extend the previous calculations in |7( 
by deriving also the so-called "non-derivative" contributions to the spin-dependent cross 
section. We have found that these combine with the "derivative" pieces into a remarkably 
simple structure. 

Using our derived cross section, we have also made first phenomenological studies, using 
the E704 fixed-target and the latest preliminary RHIC (STAR and BRAHMS) data. We 
have found that a simultaneous description of all these data is possible, albeit at a more 
qualitative, than quantitative, level, with the RHIC data overall better described. The 
"non-derivative" contributions we have calculated are of moderate importance. 

We have finally made predictions for a number of other single-spin observables at RHIC, 
in particular for the ^-dependence of for 7r° production, for scattering at 62.4 GeV, and 
for the asymmetries for single-jet and prompt-photon final states. 

For the future, it will be desirable to extend our work in a number of ways. Regarding the 
theoretical framework, one should eventually include also purely gluonic higher-twist corre- 
lation functions. These are expected to be of particular relevance for the spin asymmetry 
at mid-rapidity, which was found experimentally to be small |3]. Also, we have so far only 
considered the "soft-gluon" contributions to the spin asymmetry, for which the gluon in the 
twist-three quark-gluon correlation function is soft. As we mentioned earlier, there are also 
in general "soft-fermion" contributions. These involve, among other things, the functions 
To,,f(x, 0), rather than the T ai p{x, x) that we found for the soft-gluon case. The soft-fermion 
contributions have their own hard-scattering functions and may make a significant contri- 
bution to the spin asymmetry as well. Further points of interest will be the evolution of 
the functions T aF (x,x), and the detailed study of similarities and differences between our 

spin asymmetry in the process p^p — > nnX 
has been considered in the context of gauge links in hard-scattering processes. This could 
be achieved for example by a study of the single-spin asymmetry for two-pion or two-jet 
production in the framework of [2j that we have used here. 

Regarding phenomenology, our studies so far have a more illustrative character. With 
new experimental information arriving from RHIC, however, we will be entering an era 
where detailed global analyses of the data on will become possible. 
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Appendix A: Hard-scattering functions for inclusive-hadron production 



In this Appendix we list the hard-scattering functions relevant for single-inclusive hadron 
production. For each partonic channel, we give the functions H^ b ^ c and H^ b ^ c , which are 
to be used in Eq. ([28)1 . We also present the corresponding unpolarized cross sections H^ b _^ c 
for Eq. flU. We have: 
qg — > qg scattering: 
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19 ~~ > 91 scattering: 
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(36) 



(37) 



qq' — > qq' scattering: 



qq' — > q'q scattering: 
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qq — > qq scattering: 
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where C F = {N 2 - l)/2N c = 4/3. 
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Appendix B: Hard-scattering functions for direct-photon production 



In this Appendix we list the hard-scattering functions relevant for single-inclusive prompt- 
photon production. In this case, there are only initial-state contributions H^ b ^ in Eq. ()28jl. 
For convenience, we also again give the unpolarized contributions H\ 
19 ~ > 11 scattering: 

„2 



lab-^r We nave: 
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(43) 



11 ~ ¥ 19 scattering: 



, 2C p 
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